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^ ■ On a periodic planar graph whose edge weights satisfy a certain simple geometric 

, condition, the discrete Laplacian and d operators have the property that their deter- 

t^H \ minants and inverses only depend on the local geometry of the graph. 

' We obtain explicit expressions for the logarithms of the (normalized) determinants, 

as well as the inverses of these operators. 

We relate the logarithm of the determinants to the volume plus mean curvature of 
J> ■ an associated hyperbolic ideal polyhedron. 

In the associated dimer and spanning tree models, for which the determinants of d 
(•~^ \ and the Laplacian respectively play the role of the partition function, this allows us to 

CN ' compute the entropy and correlations in terms of the local geometry. 

■ In addition, we define a continuous family of special discrete analytic functions, 
^) ', which, via convolutions gives a general process for constructing discrete analytic func- 

' tions and discrete harmonic functions on critical planar graphs. 

■ 1 Introduction 

Let ^ be a graph and c a nonnegative weight function on the edges of Q defining the 
edge conductances. We define the Laplacian operator A : C'^' C'^l by 

A/(n) = ^ c{uw)U[u) - f{w)) 

where the sum is over vertices w adjacent to u. This is one of the most basic and useful 
operators on Q. 

On the superposition of G and its planar dual G* , we can define another operator, 
the d operator, which defines "discrete analytic functions" (those which satisfy df = 0). 
In fact we will define d on more general bipartite planar graphs, see section ^ below. 
In the form we use here the definition comes from a statistical mechanics model, the 
dimer model: the determinant of the d operator is the partition function of the dimer 



X 



model (see section pj)]) . 

We show that under a certain geometric condition, called isoradiality, on the em- 
bedding of a planar graph, the local structure of these two operators A and d takes on a 
surprising simplicity. In particular one can write down explicitly the determinants and 
inverses of these operators, which only depend on local quantities: the determinants 
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only depend on the edge weights (not the combinatorics) and the value of the inverse 
on two vertices vi,V2 only depends on a path between the vertices. 

In particular for the determinants we have the following results. For the definitions 
see below. 

Theorem 1.1 Let Qt be an isoradial embedding of a periodic planar graph and c the 
associated weight function (giving the conductances on the edges). Then the logarithm 
of the normalized determinant of the Laplacian is 

logdetiA = i V i(L(0) + L(^-0))+-logtan0, 
IN ^ — ' vr \ z / vr 

edges e 

where N is the number of vertices in a fundamental domain, the sum is over the set of 
edges in a fundamental domain, c(e) = tan(0(e)) and L is the Lobachevsky function, 

L{x) = - [ log2sintdt. (1) 



Jo 

The determinant of the Laplacian is useful since it computes the partition function 
for several other statistical mechancial models, among them the spanning tree model 
(see e.g. Q) and the discrete Gaussian free field. For the d operator we have 

Theorem 1.2 Let be an isoradial embedding of a periodic bipartite planar graph 
and V the associated weight function on edges. Then the normalized determinant of the 
discrete d operator satisfies 

logdeti9=i V -L(0) + -log2sin0, 
JN ^-^ IT vr 

edges e 

where v{e) = 2 sin 0(e), is the number of vertices in a fundamental domain, and the 
sum is over the edges in a fundamental domain. 

For results on and explicit expressions for the inverses of these operators, see The- 



orems 4.2 and 7.1 below. 

Similar statements can be made for nonperiodic graphs, on condition that the 
right-hand sides makes sense; for this one needs a unique ergodicity assumption on the 
translation action on the graph. We will restrict ourselves to periodic cases here. 



To a graph with critical weight function v as in Theorem 1.2 is naturally as- 
sociated a 3-dimensional ideal hyperbolic polyhedron (with infinitely many vertices): 
it is the polyhedron P in the upper-half space model of 'M? whose vertices are those 
of Q*j^, the dual graph oi Qd. We give a geometric realization of a discrete analytic 
function as a deformation of the associated embedding of the graph, alternatively as a 
deformation of P. Furthermore there is a relation between the above determinants and 
the volume of the polyhedra, which identifies the entropy of the corresponding dimer 
model to the volume, and the mean energy to the mean curvature. 

These results were inspired by the work of and discussions with Christian Mercat 
who used isoradial embeddings ("critical" embeddings in his terminology) in the 
study of discrete analytic functions and the Ising model. Isoradial embeddings appar- 
ently first appeared in the work of Duffin [^. Our definition of "elementary" harmonic 
functions, see section |9|, also appeared independently in [0], where a more detailed 
study of these functions is undertaken. 

The outline of the paper is as follows. 
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Section |2|. We define "isoradial" embeddings and critical weight functions for 
planar graphs. 

Section |^. We define the B operator, discrete analytic functions, and show in what 
sense they give perturbations of the embedding. We also prove the existence, 
unicity and asymptotics of the inverse d~^, as well as giving an explicit local 
formula. 

Section |5[ We define the normalized determinant deti B and prove Theorem |1.2| . 

Section |^. We point out the relation between the Laplacian and our B operator 
in the case of graphs Qt- 

Section |^ We compute the Green's function on Qt and its asymptotics. 

Section |^. We prove that the set of isoradial embeddings of a graph is convex (and 
usually infinite dimensional), and logdeti A and logdeti B are concave functions 
on it. 

Section |9|. We discuss a general construction of discrete analytic and discrete 
harmonic functions on critical graphs. 



Section 10. We discuss connections between the d operator and the dimer model. 
Here we relate the volume and mean curvature of the polyhedron P with the 
entropy and mean energy of the dimer model. 



2 Critical weight functions 
2.1 Polyhedral embeddings 

Let ^ be a planar graph, and Q* its planar dual. For a vertex v € Q we denote v* its 
dual face of Q*. Similarly for an edge e we denote by e* its dual edge. 

A polyhedral embedding of G is an embedding of Q in the plane such that every 
face is a cyclic polygon, that is, is inscribable in a circle. A regular polyhedral 
embedding is a polyhedral embedding in which the circumcenters of the faces are 
contained in the closures of the faces. 

A polyhedral embedding of a graph Q defines a map from Q* , the dual of G, to the 
plane, by sending a dual vertex to the circumcenter of the corresponding face of G, 
and a dual edge to the segment joining its endpoints. If the embedding of G is regular, 
this map will be an embedding of G* unless some circumcenters are on the boundary 
of their faces, in which case some dual vertices may have the same image in the plane. 
In other words some edges of G* may have length 0. 

The term polyhedral refers to the fact that a polyhedral embedding gives rise to 
a three-dimensional ideal hyperbolic polyhedron, in the following way. The upper 
half-space model of hyperbolic space is identified with {{x,y, z) ^ M.^\z > 0}. A 
polyhedral embedding of G on the xy-plane is the vertical projection of the edges of an 
(infinite) ideal polyhedron P{G) whose vertices are the vertices of G- If the polyhedral 
embedding is regular then P{G) is convex; the dihedral angle at an edge is the angle of 
intersection of the circumcircles of the two adjacent faces. This polyhedron will play a 
role in section 
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2.2 Isoradial embeddings 



An isoradial embedding is a regular polyhedral embedding in which all circumcircles 
have the same radius. 

If the embedding of G is isoradial then the embedding of G* will also be isoradial, 
with the same radius. We will always take the common radius to be 1. To each edge e 
of Q we associate a rhombus R(e) whose vertices are the vertices of e and the vertices 
of its dual edge. The rhombus is therefore of edge length 1 and the angle at a vertex 
of Q is in [0, vr]. See Figure |l|. 



Figure 1: Isoradial embedding of Q (thick lines) and associated rhombi (thin lines). 

A chain of rhombi is an bi- infinite sequence of rhombi {. . . , Xj, Xj+i, . . . } with Xj+i 
adjacent to Xi along the edge opposite the edge along which Xi-i and Xi are adjacent. 
The rhombi in a chain all contain an edge parallel to a given direction. When Q is 
bipartite, the vertices of ^ in a rhombus chain are alternately black and white, so that 
a chain has a "black" side and a "white" side. 

Note that a rhombus chain cannot cross itself (since it is monotone) and two rhom- 
bus chains cannot cross each other more than once, since at the crossing point the 
rhombus which both chains have in common has orientation determined by the rela- 
tive orientations of the two chains. 

Which planar graphs have isoradial embeddings? Here is an explicit conjecture. A 
zig-zag path in a planar graph is an edge path which alternately turns right and left, 
that is, upon arriving at a vertex it takes the first edge to the left, then upon arriving 
at the next vertex takes the first edge to the right, and so on. If the graph has an 
isoradial embedding a zig-zag path corresponds to a chain of rhombi. Consequently 
two necessary conditions for the existence of an isoradial embedding are that a zig-zag 
path cannot cross itself, and two zig-zag paths cannot cross each other more than once. 
We conjecture that these two conditions are sufficient as well. 

Conjecture 1 A planar graph has an isoradial embedding if and only if no zig-zag 
path crosses itself, and no two zig-zag paths cross each other more than once. 
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2.3 Critical weights 



Suppose we have a isoradial embedding of G. For each edge e of ^ define i'{e) = 2 sin 6* 
and c(e) = tan^, where 26 is the angle of the rhombus R{e) at the vertices it has in 
common with e. We refer to 9 as the half-angle of R. The function v is the weight 
function for the d matrix and c is the weight function (conductance) for the Laplacian 
matrix. A weight function or c which arises in this way is called critical. Note that 
u{e) is the length of e*, the dual edge of e. 

Note that a critical weight function determines the embedding of G and Q* . 

It is useful to allow the rhombus half-angles to be in [0,7r/2], that is, to allow 
degenerate rhombi with half-angles or tt/2. In this case the embeddings of G or G* 
can be degenerate in the sense that two or more vertices may map to the same point. 

If the half-angles of the rhombi in a chain are in [0, 7r/2] but bounded away from 
and 7r/2 we can deform the angles of all the rhombi in the chain, simply by changing 
slightly the direction of the common parallel edge. The new rhombus angles define a 
new critical weight function. Moreover one can show that every critical weight function 
can be obtained from any other by deforming in this way along chains: see section B. 



2.4 Periodic graphs 

For the rest of the paper we suppose that G is an infinite periodic graph, that is, 
the embedding (and therefore the weight function) is invariant by translates by a two- 
dimensional lattice A in M?, and the quotient t//A is finite. We denote this quotient 
Gi. 

We use periodicity in a fundamental way in defining and proving the formulas for 



the determinants. Theorems 1.1 and 1.2. For the inverses of the operators, however, 
periodicity is not necessary. As long as the graph has a reasonable regularity, for 
example if there are only a finite number of rhombus angles, then the results on the 
inverses extend (with essentially no modifications to the proofs). 



3 The d operator 
3.1 Definition 

Let Gd be a bipartite planar graph, that is, the vertices Gd can be divided into two 
subsets B UW and vertices in B (the black vertices) are only adjacent to vertices in 
W (the white vertices) and vice versa. Suppose that Gd has an isoradial embedding. 
Let be a critical weight function on Gd- 

We define a symmetric matrix 5 indexed by the vertices Gd as follows. If vi and 
V2 are not adjacent d{vi,V2) = 0. U w and b are adjacent vertices, w being white and 
b black, then d{w, b) = d{b, w) is the complex number of length ^{wb) and direction 
pointing from w to b. If w and b have the same image in the plane then the \d{w, b)\ =2, 
and the direction of d{w, b) is that which is perpendicular to the corresponding dual 
edge (which has nonzero length by definition), and has sign determined by the local 
orientation. Another useful way to say this is as follows. If rhombus R{wb) has vertices 
w, X, b, y in counterclockwise (cclw) order then d{w, b) is i times the complex vector 
x-y. 



5 



The matrix d defines the d operator. The d matrix is also called a Kasteleyn 
matrix for the underlying dimer model, see section 10, Since B maps C'^ to 



and C'^ to C^, it really consists of two operators Obw- — > and its transpose 
BwB- C^^C^. 

As an example, note that when Qd = with edge lengths r = \/2 the d operator 
has a more recognizable form, 

df{v) := V2 (f{v + r) — f{v — r) + if{v + ri) — if{v — ri)) . 



3.2 Discrete analytic functions and perturbations 

A function f : B ^ C satisfying df = 0, that is YlbeB d{w, b)f{b) = for all w, is called 
a discrete analytic function. This generalizes the definition given in . Similarly 
a function / : — > C satisfying df = is also called a discrete analytic function. So 
a critical bipartite planar graph supports two independent families of discrete analytic 
functions. 

The simplest example of a discrete analytic function is a constant function (as the 
reader may check). 

A discrete analytic function F : B ^ C defines a perturbation of the embedding of 
as follows (this idea generalizes an example in [Q]). The perturbation is a complex 
homothety (rotation, scaling and translation) on each black face but deforms the white 
faces in a more general way. See Figure 




Figure 2: A discrete analytic function is a perturbation of the embedding oi Qd. This 
example is a special case where Qd is a. superposition of an isoradial graph and its dual, see 
section ^. 
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The perturbation is defined as follows. Fix a small e > 0. To each face b* of G'^ 
we apply the map z ^ z{l + €F(b)) + tb, where tb is a certain translation. The th are 
chosen so that if two black faces meet at a vertex, their images still meet at the (image 
of the) same vertex. This is of course only possible if for each white face the polygonal 
path formed by the translates of the images of its edges is still a closed polygon: let w 
be a white vertex of Gd, and bi, . . . ,bm its neighboring black vertices. Let e^, . . . , 
be the edges running cclw around the face w* of Q'^, so that e* is the dual edge to 
wbj. Because the edges and dual edges are perpendicular, as a complex vector we have 
e* = id{w, bj). Then dF = implies that 

d{w, bi)Fibi) + ■■■ + Biw, b^)F{b^) = 0, 

or 

e{F{bi) + ■■■ + el,F{bm) = (2) 

and since + • • • + = 0, we have 

e{{l + eF{bi)) + ■■■ + e*^{l + eF{bm)) = 0. (3) 

In conclusion we can find so that each white face is again a closed polygon. 

The inverse d~^{w,b) which we define in the next section is a discrete analytic 
function of b for each fixed w, except for a single singularity: it satisfies 

J2diw',b)d-\w,b) =5^'H = 
fees 

Therefore for a fixed w it defines a perturbation of the embedding except that the 
polygon w* is not closed: at the face w* the right-hand side of (||) is 1 not 0. So we 
can therefore define only a map from the branched cover of Q"^, branched around w, to 
C; it has a non-trivial holonomy around w. See Figure ^ for d^^ in the case Gd = 

4 The operator 

For critical bipartite graphs we define to be the unique function satisfying 

1. aa-i = Id 

2. d^^{'w, b) ^ when |5 — — > oo. 

We first show unicity, then existence in the following section. 

4.1 Unicity 

Theorem 4.1 There is at most one satisfying the above two properties. 

Proof: We first show that for all w,b the argument of d~^{w,b) is well defined up 
to additive multiples of vr. To do this we conjugate 5 by a diagonal unitary matrix to 
make it real. This conjugation simply involves multiplying the edge weights on edges 
incident to each vertex by a constant (a constant depending on the vertex). 

Fix a vertex vq of Gd- Let v be any other vertex of Gd and vo,vi, . . . , V2k = v he 
a path of rhombus edges from vq to v. Each edge vjVj^i has exactly one vertex of 



f iiw^w' 
1 1 ii w = w' 
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Figure 3: The function d ^{w, ■) on the graph 7? defines a perturbation of the standard 
embedding (here we took the e of to be 3.5). 
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Ge> (the other is a vertex of Gf))- Direct the edge away from this vertex if it is white, 
and towards this vertex if it is black. Let e'"^ be the corresponding vector (which may 
point contrary to the direction of the path from vq to v). We multiply all the edges 
incident to v with weight Sy = ibe~2^"^ (either choice of sign will do). Note that 
these multipliers are well-defined up to ±1: the product of the multipliers on a path 
surrounding a rhombus is ±1. Let S be the diagonal matrix of these weights (with an 
arbitrary choice of signs). Then we claim that each entry of S*dS is real: if w,b are 
adjacent in Gd, and e*^,e'"'^ are the edge vectors of the rhombus edge R{wb), so that 
d{w,b) =i{-e^^ + e'^), then 

S*dS{w,h) = slsbd{w,b) = ±4(s^e-^(^+<^)/2)(^(_gie ^ 

We now claim that S*dS has at most one inverse tending to zero at infinity. If 
Fi,F2 are two such inverses, then Fi — F2 satisfies S*dS{Fi — F2) = 0. Moreover by 
reality of S*dS we can suppose Fi — F2 is real (otherwise take the real or imaginary 
parts separately). 



Let wq be a fixed white vertex of Gd- As in section 3.2 we associate to each black 
vertex b a polygon similar to the face b* of Gj^- That is we associate a translate of 
[S{Fi — F2)S*{'Wo, b)] ■ b* . These translates are chosen so that if polygons b\, &2 meet at 
a vertex then their images meet at the the image of the same vertex. For each white 
face w* of Gf), the images of its edges all have the same direction: if 6 is a neighbor of 
w then the dual edge (wb)* = e*^ — e**^ is mapped to the vector 

(e^e _ ^^<^>•^s{F^ - F2)S*{b,wo) = {e^' - e'^){s^e-'^'^'^^'^)sl^, 

which is sJoqSm, times a real number. Hence up to ±1 the direction is independent of 
the choice of neighbor b oi w. The image of these edges is therefore a segment. Thus 
we have a map from G*£) to the plane, such that each black face is mapped to a similar 
copy of itself (positively oriented) and each white face is mapped to a segment. Each 
such segment has a neighborhood in of its interior covered by the images of its 
neighboring black faces. If {Fi — F2){b,W()) is non-zero at some 6, the image of this 
map is open, since the image extends across the boundary of each black polygon. But 
this is impossible if Fi — F2 tends to zero at 00. □ 



4.2 Existence 

For a critical weight function on a bipartite graph Gd have the following explicit 
description of d~^. 

Let Wq be a white vertex. Define for every vertex v a rational function fv{z) as 
follows. Let Wq = vo,vi,V2, ■ ■ ■ ,Vk = V he a path in the rhombus tiling from wq to v. 
Each edge VjVjj^i has exactly one vertex of Gd (the other is a vertex of ^|)). Direct the 
edge away from this vertex if it is white, and towards this vertex if it is black. Let e*"^ 
be the corresponding vector (which may point contrary to the direction of the path). 
We define inductively along the path, starting from 

U = 1- (4) 
If the edge VjVj^i leads away from a white vertex or towards a black vertex, then 
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else if it leads towards a white or away from a black, then 

A+i=A-(^-e'"0- (6) 

To see that fv{z) is well-defined, it suffices to show that the multipliers for a path 
around a rhombus come out to 1. For a rhombus R{wh) with vertices w,x,b,y such 
that the vector wx is e*^ and the vector wy is e**^, we have 

Uz) = fy{z){z - e^'t') = h{z){z - e%{z - e^) = Uz){z - e^'). 

This shows that / is well-defined. 

Any function fv{z) satisfying ®,(^ (but not necessarily (§)) is a discrete analytic 
function for Q^i, see Theorem below. We call such functions special discrete 
analytic functions. 

For a black vertex h the coupling function d~'^{wQ, h) is defined roughly as the sum 
over the poles of fb{z) of the residue of fb times the angle of z at the pole. However 
there is an ambiguity in the choice of angle, which is only defined up to a multiple of 
27r. 

To make this definition precise, we have to assign angles (in R, not in [0,27r)) to 
the poles of fb{z). We work on the branched cover of the plane, branched over wq, so 
that for each black vertex b in this cover we can assign a real angle to the complex 
vector b — wq, which increases by 27r when b winds once around wq. Each pole of fb 
corresponds to a rhombus chain separating b from wq, with b being on the black side 
of the chain. Because the chain is monotone we can assign an angle to the common 
parallel (recall that this points from the white side to the black side of the chain) which 
is in {9q — tt, 6*0 + tt). This assigns a real angle to each pole of fb- 



Theorem 4.2 We have 



d-\wo,b) = -^ e-Res,^,.s{fb), (7) 



2Tri 

poles e'* 

where the angles G M are chosen as above for some lift of b. This can be written 

4^ /^A(^)log^rf^, 

where C is a closed contour surrounding cclw the part of the circle {e*^ | ^ G [^o ~ + 
€,9o + TT — e]} which contains all the poles of fb, and with the origin in its exterior. 

Proof: Let F{b) denote the right-hand side of (^). Note that if we choose a different 
lift of b, the angles all change by a constant multiple of 2tt; however the sum of all 
residues of / is zero since fbiz) has a double zero at oo. Therefore the right-hand side 
is independent of the lift of b. 

We need to show that X^^g^ 5(t(;, 6)F(6) = 6wg{w), and F{b) tends to zero when 
b ^ oo (see section 4.1). 



Take a white vertex w, w ^ wq, and let bi,...,bk be the neighbors of w. Let 
gidj ^^i4>j )-,g ^j^g edges of R{w,bj), so that d{w,bj) = i(— e*^J + e*'^j). For each pole e*^ 
of some fb , the value of the angle 9 is the same for each of the bj (they are part of the 
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same chain for which e*^ is the common parallel). Let Cq be a small loop around e*^ 
in C. We have 

k 

Y,d{w,b)F{b) = Y,i{-e'^^ +e''^^)F{bj), 

b&B j=l 

k 

i=l poles e^* 

k 



poles e^e i=l 



27r 27ri 

poles e'* 

However when w = wq we have 



'27T ^ 2i:i Co ^Az- e^'^){z - e^^^) 

poles J=l ^ ^ 
_^ Q r / V— -V 1 1 

'2^ E 2^ '^'"^^^ ^ (z-e^<^0 ~ U-e^^O 

poles " \j=l ^ ' ^ ' 

- E 



dz 



Odz = 0. 

Ce 



Fit,) = - J- ( ^-%) 



so that ^ ^ 

E z(-e^^^ + e^^^ ) = l-j^O, - c^, = 1, 

since the angles increase by 27r around wq- 

Finally to show that F{b) tends to zero see Theorem |4.3| below. □ 



4.3 Asymptotics of 5 ^ 

Let w = vo,vi, . . . , V2k = & be a path of rhombus edges from w tob which crosses each 
rhombus chain at most once. Such a path exists since two chains cross each other at 
most once. Every vertex V2j is a vertex oiQo^ and the edges coming to and leaving V2j 
contribute a pole and a zero to /{, (we are assuming fw{z) = 1). So fb{z) has the form 

■^^^ ' [z - e^^i)(z - e^^2) JLL _ gift) • 

By periodicity of Qd-, the angles /3j are all in [6q — -k + e,9Q + iT — e] for some fixed e > 
independent of b. 

Theorem 4.3 We have 



2ii \b — w {b — w) J 2tt \{b — w)^ ^{b — w)^ J \\b — ■w\'^ 
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where 7 = /b(0) 
above notation. 



-iiei+92) 



i{-f3j+aj) 



and ^2 



+ e^^^^+Y.e 



Proof: By Theorem 4.2, d ^{w,b) is given by an integral 



(27r)2i 



fb{z) log zdz 



c 



id 



on a simple closed curve C which surrounds in a cclw sense the set {e 
[^0 — vr + e, ^0 + ^ ~ e]} and has the origin in its exterior. In fact since / tends to zero 
at 00 we can homotope C to the curve from 00 to the origin and back to cxd along the 
two sides of the ray R from the origin in direction + t^- On the two sides of this ray, 
logz differs by 27ri: it is 27rz less on the cclw side, and the integral therefore becomes 



27r 



fbiz)dz, 



R 



the integral being from cxd to along R. 

We show that when |6 — w| is large the only contributions to this integral come from 
a neighborhood of the origin and from a neighborhood of 00. Suppose without loss of 
generality that Oq = 0, so that R is the negative real axis. Define 



fe-i 



i/3j iaj 



N =\b-w\ = |e*^i + e'^^ + ^ e'^^ - e'"^ I 



For small t we have 



-if3 



giving 



A(t) = 7exp((6-u;)t + e2^ + 0(t3; 



(8) 



Since f3 G [— vr + e, tt — e], we have cos /5 > 1 — 5 for some fixed 6 = 6{e). For all 
i < we can use the bound 



which gives 



t 



oif3 



< 



< 



't2 + 1 - 2tcosa 
t2 + 1 - 2t cos (3 



'1 + 



exp 



2t(cos (3 — cos a) 
t2 + 1 - 2t(l - 6) 

t(cos (3 — cos q) 
t^ + 1- 2t{l - 6) J ' 



\h{t)\ < 



|(t-e^^i)(t-e^«2)| 



< gt(JV-2)/(l+t2-2t(l-5))^ 



(9) 
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For large t we have 



which gives 

= e(''— )/*+«2/(2t')+0(*-'). (10) 

When —I/^/N < t < 0, we can use the approximation (|8|), when — A^^/^ < t < 
1/y/N we can use (P) and when t < —N^/"^ we can use (|lO|). 
The small-|t| integral gives 

The intermediate integral is negligible, and the large-t integral gives 



□ 



5 Determinants 

For a finite graph Q we define the normalized determinant deti M of an operator 
M: ^ to be 

detiM-'^^'ldetMlViei, 

where |^| is the number of vertices of Q. 

For an operator M on a finite graph, det M is a function of the matrix entries 
M{i,j) which is linear in each entry separately. In particular for an edge e — I J , as a 
function of the matrix entry M{i,j) we have detM = a + PM{i, j), where /? is (— 1)*^-' 
times the determinant of the minor obtained by removing row i and column j. That 
is, 

a(det M) 
d{M{i,j)) 

or 



M-^(j,i) - detM, 



a(logdetM) 



dM{i,j) 



M-'{j,i). (11) 



Suppose that Q is periodic under translates by a lattice A. Let Qn = G/nA, the 
finite graph which is the quotient of G by nA. Now if we sum (^) for all A-translates 
of edge ij in Gn, and then divide both sides by \Gn\, it yields 

0(logdet,M)^ 1 

dM{i,j) \Gi\ ^' ^ ' 

where M{i,j) is now the common weight of all translates of edge ij, that is, the left- 
hand side is the change in logdeti M when the weight of all translates of ij changes. 
Note that this equation is independent of n. 
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We can now use (^) to define the normalized determinant of d on an infinite 
periodic graph in terms of d~^. Rather, it defines the change in the normahzed 
determinant. To complete the definition we need to define the normalized determinant 
of the 3 operator at some initial condition. If the half-angles of a graph degenerate by 
tending to or 7r/2, the edge weights tend to either or 2. The edges with weight 
contribute nothing, the edge of weight 2 each are defined to contribute \/2 to the 
normalized determinant. Thus we define the limit of logdeti d to be log 2, where 
n is the number of edges per fundamental domain with weight 2. This makes sense 
since under this degeneration, removing the zero weight edges, the graph consists of 
copies of Z, each consisting of weight 2 edges. 



Proof of Theorem 1.2; Suppose Qd is periodic with lattice A. Take a rhombus 
chain and all its A-translates in Qjy. Let (po be the smallest half-angle of a rhombus on 
these chains, and (pi the largest half-angle: they exist by periodicity. If (po < Tr/2 — (pi, 
deform these chains by changing the common parallels to decrease (po to zero: at that 
point the corresponding edge weight is zero and we can remove the edge completely 
from Qo, simplifying the graph in the process. If on the other hand 7^/2 — (pi < (pQ, 
deform the chain by increasing (pi to 7r/2. At this point the rhombus Ri of angle (pi 
has flattened out to become a segment. We now repeat with another chain, but in 
future deformations we always keep Ri as a segment: the angle (pi remains tt/2. We 
need to show that this is always possible. If we deform along a chain or union of chains 
whose common parallels do not contain an edge of Ri, the deformation will not affect 
Ri. If an edge of Ri is in the common parallel of a chain being deformed, we deform 
simultaneously and in the same direction both chains which pass through Ri. So a 
deformation involves changing the common angle of parallels of chains, all of which 
have a parallel in the same direction. We can deform in this way because for any set 
of chains having parallels in the same direction, there is a deformation which tilts all 
the common parallels and is a translation on each complementary component of these 
chains. We continue deforming until all angles are zero or 7r/2. We have proved that 
any rhombus tiling of the plane can be deformed through rhombus tilings to flatten out 
all the rhombi: the rhombus vertices all map to Z, the rhombus edges map to edges of 
Z, and the edges oi Qd either have length 2 (and weight zero) or length (and weight 
2). 



We can compute the change in logdeti 3 as in (12) under such a deformation since 
it only depends on the angles of the rhombi in the chains. 

To complete the calculation it suffices to compute the contribution for a deformation 
along a single chain. In fact we can treat each rhombus of the chain separately. Suppose 
we deform a chain whose common parallel is e^*^. Suppose the rhombus R{wb) has 
edges e*" and e*("+^^) leading away from w. Then the term in 

dlog deti B 



de 



which involves this rhombus is 



rf(logdeti9) didiw,b)) ^ .1 j^. didiw,b)) 

' d{d{w,b)) ' de ' de 



L ( '^1 ^ i.fif_ei(°+2e) , 

2m Ve^("+2e) -e^°; de^ ^ 

TTsin^ 
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Suppose in the above deformations this rhombus has half-angle going to zero. Inte- 
grating the real part of this for 9 from to </> yields 

1 (j) 1 

— 9 cot 9d9 = — log sin (p / log sin 9d9 

Jo TT Jq 



(integrating by parts), which is 



- log 2 sin (/) + -L{(f)). 
TT vr 



On the other hand if the rhombus half-angle goes to vr/2, we should integrate from 
vr down to (j). This yields 

log 2 + - log 2 sin (/) + -L((/)), 
2 vr vr 

since when (p = 7r/2 we have 

— log 2 sin (/) + — = log 2. 
vr vr 2 

This completes the proof. □ 

Note that for in its standard embedding this gives the classical result of Kasteleyn 
m, except for an extra factor ^ log 2 per edge since our edge weights are i^(e) = \/2, 
not 1. 

A more classical way of defining the normalized determinant is via an exhaustion 
Gd by finite graphs. On an infinite graph Q^o which is a limit of a growing sequence 
of finite graphs Goo = ^i^n-fcx>Gn, and supposing that an operator M^o on Goo is a 

de f 

limit in an appropriate sense of Mn on Gn, one would wish to define deti Moo = 
lim„^oo deti M„, assuming this limit exists. Unfortunately it is not so easy to prove 
that the limit exists in the case of the operator d, even when the graph is with 
periodic weights: see M. That is why we chose the previous definition instead. 



6 Relating the Laplacian to d 

Let Gt be a planar graph (not necessarily bipartite) with an isoradial embedding in the 
plane. To an edge e with rhombus angle 29 £ [0,vr] we associate a weight c(e) = tan^. 
This weight c(e) is the conductance in the underlying electrical network. This is the 



situation considered in [|10|. 



Given Gt we define another graph Gd to be the superposition of Gt and its dual, 
that is, Gd has a vertex for each vertex, edge and face of Gt, and an edge for each 
half-edge of Gt and each half dual-edge as well. See Figure The graph Gd '^s a 
bipartite planar graph, with black vertices of two types: vertices of Gt and vertices of 
Gtj and white vertices of Gd correspond to the edges of Gt- 

An isoradial embedding of Gt gives rise to an isoradial embedding of Gd- each 
rhombus of Gt with half-angle 9 is divided into four congruent rhombi in Gd, two of 
half-angle 9 and two of half-angle ^ — 9: see Figure ^. 

An edge of Gt of rhombus half-angle 9 has weight c(e) = tan^. This defines weight 
tan^ on the two "halves" of each edge in Gt U G^ (which are edges of Gd), and the 
dual edges are naturally weighted 1 (as in [^). 
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Figure 4: Isoradial embedding of Qt (black vertices, thick lines) and its dual (grey 
vertices, thin lines), and associated rhombi of Qd (dotted lines). 

The B operator on Qd and the Laplacian on Qt are related as follows. Let D be 
obtained from B by multiplying edge weights oi Qd around each white vertex (coming 
from an edge of Gt of half-angle 6*) by l/(2\/sin0cos0), that is, D{w, b) = Sd{w, b)S* 
where S is the diagonal matrix defined by S{w,w) = l/{2^/sin9{w) cos9{w)) and 
S{b,b) = l. 

Then an edge of coming from a "primal" edge of Qt has weight VtanO for D, 
and an edge coming from a dual edge has weight l/VtsmO for D. (These edges have 
weight 2 sin ^,2 cos respectively for d.) 

Lemma 6.1 Restricted to vertices of Gt we have D*D = Ag^,. Restricted to faces of 
Gt, we have D*D = Ag* , the Laplacian on Gt- 

Proof: Let b & Gt have neighbors 5i, . . . , 6^ and let 6i be the half-angle of edge bbi. 
Then in Qd, b has neighbors wi, . . .w^ with these same half-angles, and D{b,Wj) = 
aj tan Qj where \oLj \ = 1. 
We have 

k k 

D*D{b,b) = \aj^/taJ[e'j\'^ = J^^tan^^- = A{b,b). 

j=i j=i 
If b, b' are adjacent in Gt we have 

D*D{b,b') = {aVtane){-aVtane) = -tanO = A{b,b'). 

If 6,6' are at distance 2 vuGd but correspond to a vertex and a face of Gt-, we have 
two contributions to D*D{b,b') which cancel: in Gt let ei,e2 be the two edges from 6 
bounding the face 6', with ei on the right. Then 

D*D{b,b') = D{b,wi)D*{wi,b')+D{b,W2)D*{w2,b') 

. I -— , , — ioT I — N 

= (aiVtan6li)(^===) -f- (a2Vtan6'2)(^===) 
V tan B\ V tan Hi 

= 0. 
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Finally, if b' is at distance at least 2 in Qt-, D*D(b, h') = 0. 

The proof for is identical. □ 



7 Green's function 

As we did for we can give an explicit expression for the Green's function on a 
critical graph. 

The Green's function G{v,w) is defined to be the function satisfying 

1. AG{v,w) = 5y{w) (the Laplacian is taken with respect to the second variable), 

2. G{v, w) = 0(log(|t; - w\)), and 

3. G{v,v) = 0. 

We define for each vertex v of Gt and a rational function gv{z) of 2 as follows. 
We define gy inductively starting from g^g = ^ and if v,v' are adjacent vertices of a 
rhombus then 

9v'=gv-j—[Q, (13) 

where e*^ is the complex vector from v to v'. Now g^ is well-defined since the product 
of the multipliers around a rhombus is 1. 

Theorem 7.1 

G{vo,vi) = ^ X] ^ ■ ^^^2=e'«(5«;i(^)) = / 9viiz)logzdz, 

poles e'-^T^O ^ 



where the angles 9 and curve G are chosen as in Theorem J^.i. Forvo G Gt andvi G Gt 
this formula defines i times the harmonic conjugate of G. 



Proof: The proof is similar to the proof of Theorem O. Clearly G{v,v) = by 
definition. To show that AG{v,w) = 6y{w), we show that for each edge the discrete 
Cauchy-Riemann equations are satisfied (||lO|), that is, we show that for an edge vv' 
with half-angle 9 and with dual edge ww' , we have 

{G{vo,v') — G{vo,v))i tan 9 = G{vq, w') — G{vq,w), 

assuming ww' is obtained by turning cclw from vv'. It suffices that for each pole e*"^ 
we have 

idv'iz) - gv{z))i tan 9 = g^'{z) - g^iz). (14) 

By (^), we have 

z + e^ z + e^ (z + e^°)(z + e'/^) 

gy,,{z) = gy{z)-—-, gUz) = gy{z)-—^, g„,(z) = ^^z)^-—-^^—-^^ 

so that dH) is 

z + e'^'z + e'^ \. a-P , . ( z + e''^ z + e*^ 



9v{z) -77- -fl - 1 itan— — = gy{z) 
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This identity holds, therefore G{vo,w) is harmonic for w ^ vq. Moreover for a dual 
vertex v* adjacent to vq we have 

Qv {z) 



Z{Z - 6'-°') ' 



whose residue at e*" is 2. Thus the value of G{vq,v) is —ai/2TT. In particular for one 
cclw turn around the vertex vq, the dual function increases by — i, which implies that 
the Laplacian at vq is 1. Finally to show that G{v, w) = 0(log(|t; — w\)), see Theorem 



7.3 below. □ 
From Lemma 16. 11 there is a relation between G and D^^: 



Corollary 7.2 When edge bib'i of Qt corresponds to vertex wi oJQd, we have 

1 



Proof: By Lemma |6.1| , taking the Laplacian of D^^{wi,v) with respect to v gives 

^g^D-\wi,v) = D*DD-\wi,v) 

where /i, /2 are the two faces adjacent to ei. In particular D^^{wi, b) is harmonic as a 
function of b £ Qt except at bi and b'^, where its Laplacian is zizD*{wi, bi) respectively. 
Since D^^{wi, 62) ^ as 62 ^ C!0 this proves (p^). □ 



7.1 Asymptotics of the Green's function 

Theorem 7.3 

G{vi,V2) = --^log\v2 - vi\ -^^^^ + 0{-. 



2tt 27r \v2 — vi\ 

Note that our Laplacian is positive semidefinite, which results in the minus sign in 
this formula. 

Proof: Let vq = wqjWi, . . . ,Wk = vi be a path of rhombus edges from vq to vi. 
Let e*^J be the complex vector from wj to wj+i. As in the proof of Theorem |4.3| we 
assume without loss of generality that 6j £ [— vr + e, vr — e] for some fixed e independent 

of Vl,V2- 

Let iV = l-ui - uo|. We take r < 1/iV and > iV. We integrate 



T^i J (J z z — e*^J ' 



where C is a curve which runs cclw around the ball of radius R around the origin, from 
angle — vr to vr, then along the negative x-axis from —R to — r, then clockwise around 
the ball of radius r from angle vr to — vr, and then back along the x-axis from — r to 
-R. 
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The integral around the ball of radius r is 

^ ^ "(-l)^(l + 0{Nr)){\ogr + iB)id6 = ^"^^^^ ^"^"^ (l + 0{Nr)). 



The difference between the value of log z above and below the x-axis is 27ri, so that we 
can combine the two parts of the integral along the x-axis into the integral of 



This integral we split into a part from —R to —vN, from —\/N to —l/\/N and from 
-1/y/N to -r. 

For small t < we have 

k—l iQ . 

11 ^ _ g^idj 

j=0 

SO that the integral near the origin is 



4vr ./-i/Viv * 



2^ 

dt 



where a = 2{vi — vq), 



1 p-ar s 

(1 + 0(t7))/ -ds 



47r N J_a/VN s 

— ^1 + 0— / —ds+ ds+ — 

-l)^ _ 1 

(log(ra) + 7Euler) + 0{ — ). 



The Euler 7 constant 7Euior comes from evaluating the first two of these integrals in 
the limit ar — > and a/\fN 00. 

In the intermediate range the integral is negligible (see the proof of Theorem |4.3| ) . 
Near t = —R we have 



t - e^-^i 

so that the integral near —72 is 

dt = / ds, 



1 / /■-! - 1 , f-^ ds /-^/° 

-ds+ — + 



yj-R/a S J-R/a S 7_i 

(- log(i2/a) + 7Euler) + O(^), 
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These are precisely the same integrals as in the "near r" case, under the change of 
variable s 1/s. 

Finally the integral around the ball of radius R gives — ^ log R. The sum is there- 
fore 

-^aog« + (-1)' iog«) - (1 + (-1)')^ + o(^) 

which when k is even is 

1 1 I I 7Euler , ^/ 1 N 1 o / ^ 7Euler , r,t ^ \ 

log |.i - -ol - ^ + O(-) = --Re(.i - -o) - ^ + O(^) 

and when k is odd it is 

log + O(i-) = -f Imlog(^i - vo) + O(^). 

47r -ui — A* ZTT iV 

□ 



7.2 Determinant of Laplacian 



Using Lemma |6.1| we can compute the determinant of A in terms of the determinant 
of d. However it is as simple to compute it directly. 

We define the normalized determinant deti A to be the function of the conductances 
satisfying, for w adjacent to v, 

(?(logdetiA) 1 , ^, N ^/ ^ ^^ 2G(v,w) 
Iff J = TTTT {G{v, V - G{v, w) - G{w, v) + G{w, w)) = 

and for a graph with zero conductances, deti A := 1. 

Let Qt be a periodic planar graph, periodic under A with fundamental domain 

gT/A = gi. 

Proof of Theorem |1.1| The method of proof is the same as that of Theorem L2 . Take 



a chain of rhombi and all its translates under A, and change their common parallel, 
decreasing the smallest angle to zero. If an edge has angle which goes to zero, its 
conductance goes to as well and so it can be removed from Qt- Similarly if any 
rhombus angle increases to 7r/2, the corresponding conductance becomes infinite and 
one can remove the corresponding edge, gluing the two vertices together. In this way 
we can simplify the graph until no edges remain. 

However there is a problem which is that the determinant blows up when a con- 
ductance goes to oo. So instead we compute the change in the difference 

1 sr^ 29 

f := logdetiA — ■— -- > — log tan 

as we reduce the graph in this manner. (The quantity £ is the entropy of the underlying 
spanning tree model.) 

To compute the change under the perturbations we can treat each rhombus inde- 
pendently. For a rhombus R{vw) of angle 26, with edges from v given by e*" and 
gi(o+2e)^ a short calculation yields 

G{v,w) = 

vr tan 
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so that the part of d£/d9 involving the rhombus is l/|^i| times 

,d{Mv,w)) d [26^ \ 20 dtan9 /21ogtan^ 26 
-2G{v,w) - logtang = ^ + 



89 89 \Tr J TT tan 9 89 \ it tt tan 9 cos 

_ 2 log tan 9 

TT 

The integral from to is 



2 

TT Jo 



2 

logtan^# = -(L(^) + L(7r/2 - ^)), 

TT 



where we used L(7r/2) = 0. Similarly the integral from 7r/2 down to cf) is 

-iL{cP) + L{7r/2-4>)), 

TT 

using L(0) = L(7r/2) =0. □ 



8 Convexity 

Proposition 8.1 When parametrized by the rhombus angles the set of isoradial em- 
beddings of a graph Q is convex. 

Proof: Each rhombus has angle in [0,7r]. The set of allowed angles is the subset of 
[0, 7r]l^l defined by the linear constraints that the sum of the angles around each vertex 
be 27r. It is therefore convex. □ 



Theorem 8.2 A periodic graph which has an isoradial embedding has a unique critical 
weight function v maximizing deti ^, and a unique critical weight function c maximizing 
deti A. 

Proof: Given that the set of critical weight functions is convex, and Z := deti 8 or 
deti A is continuous as a function of the rhombus angles, we know that Z attains its 
maximum. To show unicity, it suffices to show that log Z is a strictly concave function 
of the rhombus angles. Suppose there are n rhombi in a fundamental domain. Then 
log Z is the sum over the rhombi in the fundamental domain of a function depending 
only on the rhombus angle. Moreover, the functions 

fD{9) = -L{9) + -log2sm9 

IT TT 

and 

frie) = - (m+L{^ - 9)) + - log tan ^ 

are strictly concave as a function of 9. So we can view logZ as a function on [0,7r]" 
(the set of rhombus angles), subject to the linear restrictions that the sum of the 
angles around each vertex be 27r. Thus log Z is the restriction to a linear subspace of a 
strictly concave function on [0,7r]". As such it is strictly concave and so has a unique 
maximum. □ 

What are the geometric intepretations of these maxima? 
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9 Discretization 

Recall that a special discrete analytic function is a function fb{z) satisfying (^) and 
(|^. Note that it is defined by its value at any vertex. Other discret analytic functions 
can be built up from these by convolutions: 

Theorem 9.1 Let fb{z) he a special discrete analytic function and jj, he any measure 
in the plane. Then the function F{b) = J fh{z)d^{z) is a discrete analytic function. 

Proof: The proof was already given, in the proof of Theorem we repeat it here. 
Let w he a, white vertex with neighbors hi, . . . ,bk, so that R{wbj) has edges e*^^ and 



We have 



k 

^d{w,h)F{b) = ^i{-e^^^ +e''l'^)F{b 

b(^B j=l 
k 



/ ^-^'^ - ] ^^^^^ 

/ 



M^i{z) = 0. 



n 



There is corresponding version for discrete harmonic functions, which also appear 
in a different form in [^]: 

Corollary 9.2 On Qt, let hy^lz) be given and if hj]{z) is defined and v' is adjacent to 
V, define 

{z + e'%{z + e'^) 

Ki{z) = K[z)- 

[z — e^^)[z — e^f) 

where 9 and (f) are the angles of the edges of rhombus R[vv'). This defines hy{z) for 
all V G Gt, o^nd H(v) = J hv{z)dfi{z) is discrete harmonic. 



10 The dimer model 
10.1 Definitions 

The dimer model on planar graphs was initiated by Fisher, Kasteleyn and Temperley 
||6|, |l^, who "solved" the model on 7? by finding a closed-form expression for the 
partition function. 

A dimer covering, or perfect matching, of a graph ^ is a set of edges M of ^ which 
covers every vertex exactly once, that is every vertex is an endpoint of exactly one 
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element of M. Let u: E ^ [0, oo) be a weight function on the edges of Q. To an 
edge e we associate an energy — logz/(e); this defines a natural probability measure ^, 
the Boltzman measure, on the set of all dimer coverings where the probability of 
a configuration is proportional to the product of its dimer weights (the exponent of 
minus the sum of the energies), ^{C) = rieec ^(^)- dimer model is the study of 
this measure. 

The dimer partition function is by definition 

MeM 

The partition function per site is by definition Z = i/)-*^/!^! , where \Q\ is the 
number of vertices of Q. 

For infinite graphs one can define the partition function per site by taking limits of 
the partition functions per site on finite pieces: if {Gn} is an exhaustion of Q by finite 
subgraphs, one can define 

Z= lim Z(g„,i/)i/l^"l, 

n— »oo 

when this limit exists. Usually this limit will depend on the sequence {^n}: see e.g. 
i- 

One can also define a probability measure fi on M{G) which is a weak limit fj, = 
lim^^oo fJ-n where fin is the Boltzmann measure on M{Gn), assuming this limit exists. 
This limit is even more dependent on the sequence than the partition function per site. 

One general situation where these limits are known to exist is when Gd is of the 
form Qt U as described in section |6[ For a periodic Gt the Green's functions are 
known to converge to the corresponding Green's function on the plane as defined in 
section ^, see 0. From corollary 15 we can compute the operator from the Green's 



function. In fact in this situation one can exhibit a measure-preserving isomorphism 
from the dimer model on Qd to the corresponding spanning tree model on Qt |8|. 
In particular the partitions functions are equal. 

The utility of the d operator in the dimer model is the subject of the next section. 



10.2 Dimers on finite graphs 

A subgraph Qi of Q is said to be simply connected if it consists of the set of edges 
and vertices contained on or inside a simple polygonal path in Q. 

Theorem 10.1 Let Qi he a simply connected subgraph of a graph Q (where Q has a 
critical weight function), and di the submatrix of d corresponding to Qi. The dimer 
partition function on Qi is Z\ = a/ | det d\ \ . The probability of edges \w\h\^ . . . , Wfc^fc} 
occurring in a configuration chosen with respect to the Boltzmann measure jj, is 

{X{di{w,M)\ ^ det ((ai)-iK,5,)). 

The function d^"^ is called the coupling function or inverse Kasteleyn matrix. 

Proof: For the first statement it suffices to show that B is Kasteleyn-fiat Q. 
That is, we must show that for each face of Qi (except the outer face) with vertices 
til, fi, . . . , Um,Vm in cyclic order, we have 

arg{d{ui,vi) . . . d{um,Vm)) = arg((-l)'""^(9(?;i, n2) . . . d{vm-i,Um)d{vm,ui)). 
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(This identity implies that two dimer configurations which only differ around a single 
face have the same argument in the expansion of the determinant. By any two 
configurations can be obtained from one another by such displacements.) To prove 
this identity, note that it is true if the points are regularly spaced along a 2m-gon, and 
note that it remains true if you move one point at a time. 

For the second statement, see M. □ 



10.3 Volume and mean curvature of P{G) 



Proposition 10.2 (Milnor [14]) The volume of a hyperbolic simplex in the upper 
half space model with vertices at oo, (0, 0, 1), (1, 0, 0) and (cos 9, sin 6*, 0) is L{{tt — 6)/2), 
where L is the Lobachevsky function 

Let ^ be a bipartite planar graph with an isoradial embedding, and P{G*) the 
hyperbolic polyhedron associated to its dual (vertices of P are vertices of G*). For 
each rhombus R{wb) of Q let w, x, b, y be its four vertices; we can decompose P into 
simplices, so that the part of P over R consists of two simplices with vertices w,x,y,oo 
and b, x, y, oo. Here vertices x, y are ideal whereas w and b are in the interior of 
hyperbolic three-space (at Euclidean distance 1 from the xy-plane). Each of these 
simplices has volume L{9), where 6 is the half-angle of the rhombus. Therefore in the 
expression for the determinant of B, the sum of the terms involving L is 2^ times the 
volume of P. 

The mean curvature of a compact convex polyhedron in is by definition the sum 
over the edges of (vr — 0dih)^5 where 6'dih is the dihedral angle and i is the hyperbolic 
edge length. For an ideal polyhedron, the edge lengths are infinite, so one first assigns 
a small horosphere to each vertex and then computes the edge lengths in the exterior of 
these horospheres. This definition of course depends on the radii of the horospheres. In 
the case of a P arising from an isoradial embedding of Q we can choose for a fixed small 
e all the horospheres to be Euclidean balls of diameter e; in this case the hyperbolic 
length of an edge between points vi and V2 is exactly 

21og- 



e 

where |v2 — is the Euclidean distance from vi to V2- 

For a rhombus of half-angle 9, the dihedral angle of P at the dual edge is n — 29. 
Therefore the mean curvature per vertex of P is 

^ eiog^^ = -47rloge + -^ ^ 01og2sin0 (16) 

this last equality holds since at each vertex the sum of the half-angles is vr. Finally it 
makes sense to define the normalized mean curvature to be the above formula (^) 
without the —4tt log e term. 

10.4 Geometry and Statistical mechanics 

li Qd arises as Qt U for a periodic isoradial graph Gt, we can identify deti d with 
the partition function of the associated dimer model, and as the coupling function, 



as discussed at the end of section |10.1 
Moreover we have 
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Theorem 10.3 The volume per site of P{Q*j-,) is 2tt times the entropy per site of the 
dimer model. The normalized mean curvature per site of P{G'^) is Att times the mean 
energy per site of a dimer configuration. 

Proof: In the dimer model the probabihty of occurrence of an edge e = wh \s 
i'{e)\d-^{w,h)\ = f. Recall that the energy contribution for this edge is log 2 sin ^. 
The average energy per vertex E is therefore l/|^i| times the sum over the edges in 
a fundamental domain of ^ log 2 sin 0. By the argument of Proposition 9.1 of the 
fact that the edge correlations tend to zero implies that the entropy per vertex of the 
dimer model is precisely (logdeti B) — E, or 

' ' e 

that is, l/27r times the volume per vertex of P. □ 
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